Stationary non-equilibrium properties for a heat conduction model 
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We consider a stochastic heat conduction model for solids composed by N interacting atoms. The 
system is in contact with two heat baths at different temperature Te and T r . The bulk dynamics 
conserve two quantities: the energy and the deformation between atoms. If Te T r , a heat flux 
takes place in the system. For large N , the system adopts a linear temperature profile between 
Te and T r . We establish the hydrodynamic limit for the two conserved quantities. We introduce 
the fluctuations field of the energy and of the deformation in the non-equilibrium steady state. 
As TV goes to infinity, we show that this field converges to a Gaussian field and we compute the 
limiting covariance matrix. The main contribution of the paper is the study of large deviations for 
the temperature profile in the non-equilibrium stationary state. A variational formula for the rate 
function is derived following the recent macroscopic fluctuation theory of Bertini et al. 

PACS numbers: 44.10.-K, 05.60.-k, 63.10.+a, 66.70.+f 



I. INTRODUCTION 

Understanding of the steady state of non-equilibrium 
systems is the subject of intense research. The typical 
situation is a solid in contact with two heat baths at 
different temperature. At the difference of equilibrium 
systems where the Boltzmann-Gibbs formalism provides 
an explicit description of the steady state, no equivalent 
theory is available for non-equilibrium stationary state 
(NESS). 

In the last few years, efforts have been concentrated 
on stochastic lattice gases ([!]). For these latter pre- 
cious informations on the steady state like the typi- 
cal macroscopic profile of conserved quantities and the 
form of the Gaussian fluctuations around this profile 
have been obtained ((H). Recently, Bertini, De Sole, 
Gabrielli, Jona-Lasinio and Landim proposed a defini- 
tion of non-equilibrium thermodynamic functionals via 
a macroscopic fluctuation theory (MFT) which gives for 
large diffusive systems the probability of atypical pro- 
files (0,[3) in NESS. The method relies on the theory 
of hydrodynamic limits and can be seen as an infinite- 
dimensional generalization of the Freidlin-Wentzcl the- 
ory. The approach of Bertini et al. provides a variational 
principle from which one can write the equation of the 
time evolution of the typical profile responsible of a given 
fluctuation. The resolution of this variational problem 
is, however, in general very difficult and it has only been 
carried for two models : the Symmetric Simple Exclusion 
Process (SSEP) ([3]) and the Kipnis Marchioro Presutti 
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(KMP) model (d). 

Hence, it is of extreme importance 
to identify simple models where one can test the validity 
of MFT. 

The most studied stochastic lattice gas is the Simple 
Exclusion Process. Particles perform random walks on a 
lattice but jumps to occupied sites are suppressed. Hence 
the only interaction is due to the exclusion condition. 
The only conserved quantity by the bulk dynamics is the 
number of particles. In this situation, the heat reservoirs 
are replaced by particles reservoirs which fix the density 
at the boundaries. The KMP process is a Markov process 
composed of particles on a lattice. Each particle has 
an energy and a stochastic mechanism exchange energy 
between nearest- neigbor particles (@). 

The real motivation is to extend MFT for Hamiltonian 
systems (@). Unfortunately, for these later, even the 
derivation of the typical profile of temperature adopted 
by the system in the steady state is out of range of the 
actual techniques (0|). The difficulty is to show that 
the systems behave ergodically, e.g. that the only time 
invariant measures locally absolutely continuous w.r.t. 
Lebesgue measure are, for infinitely extended spatial uni- 
form systems, of the Gibbs type. For some stochastic lat- 
tice gases it can be proven but it remains a challenging 
problem for Hamiltonian dynamics. 

We investigate here the MFT for a system of har- 
monic oscillators perturbed by a conservative noise (0, 
These stochastic perturbations are here to re- 
produce (qualitatively) the effective (deterministic) ran- 
domness coming from the Hamiltonian dynamics (PH. 
03 > 03 )■ This hybrid system can be considered as a 
first modest step in the direction of purely Hamiltonian 
systems. 

From a more technical point of view, SSEP and KMP 
are gradient systems and have only one conserved quan- 
tity. For gradient systems the microscopic current is a 
gradient ( 03 ) so that the macroscopic diffusive character 
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of the system is trivial. Dealing with non-gradient mod- 
els, we have to show that microscopically, the current is 
a gradient up to a small fluctuating term. The decompo- 
sition of the current in these two terms is known in the 
hydrodynamic limit literature as a fluctuation- dissipation 
equation (15]). In general, it is extremely difficult to 
solve such an equation. 

Our model has two conserved quantities, energy and 
deformation, and is non-gradient. But fortunately, an 
exact fluctuation-dissipation equation can be established. 
In fact we are not able to apply MFT for the two con- 
served quantities but only for the temperature field which 
is a simple, but non-linear, functional of the energy and 
deformation fields. 

The paper is organized as follows. In section[lll we de- 
fine the model. In section lHll we establish the fluctuation- 
dissipation equation and obtain hydrodynamic limits for 
the system in a diffusive scale. Section HVl is devoted to a 
physical interpretation of the fluctuating term appearing 
in the fluctuation-dissipation equation. In section [V] we 
compute the covariance of the fluctuation fields in the 
NESS by a dynamical approach and show the covariance 
for the energy presents a non-locality we retrieve in the 
large deviation functional (the quasi-potential) . The lat- 
ter is studied in section [VI] for the temperature field. 



II. 



THE MODEL 



We consider the dynamics of the open system of length 
N. Atoms are labeled by x £ {1, . . . , N}. Atom 1 and 
N are in contact with two heat reservoirs at different 
temperatures Tg and T r . Momenta of atoms are denoted 
by pi, ... ,pn and the distance between particles are de- 
noted by T"i, . . . , rjv— i- The Hamiltonian of the system is 
given by 

N 2,2 

-H N =Y, e ^ e x = ^p x = l,...,N-l 



ejv 



P% 



We consider stochastic dynamics where the probability 
density distribution on the phase space at time t, de- 
noted by P(t,p,r), evolves following the Fokker-Planck 
equation 

d t P = N 2 L*P 

Here C = A + jS + B\^ t + <SjV,T r is the generator of the 
process and C* the adjoint operator. The factor A 2 in 
front of C* is here because we have speeded up the time 
by A 2 , this corresponds to a diffusive scaling. 
A is the usual Hamiltonian vector field 



N-l 



N-l 



A =^2 (Px+i - Px)d rx + ^2 (r x - r x -i)d Pm 

x=l x=2 

+(ri - £)d Pl - (rjv_i - £)d PN 



The constant £ fix the deformation at the boundaries. 

S is the generator of the stochastic perturbation and 
7 > is a positive parameter that regulates its strength. 
The operator S acts only on momenta {p x } and generates 
a diffusion on the surface of constant kinetic energy. This 
is defined as follows. For every nearest neigbor atoms x 
and x + 1, consider the following one dimensional surface 
of constant kinetic energy e 

S* ={(Px,P*+i) £M. 2 ;pl+p 2 x+1 =e} 

The following vector field X x ,x+i is tangent to Sg 

X x . x+ i — p x +id Px ~p x d Pu:+l (1) 

so X 2 J., j generates a diffusion on (Brownian motion 
on the circle). We define 

S = 2 X x,x+i 

x=l 

Bi^Te and Bn,t t are two boundary generators of 
Langevin baths at temperature Tg and T r 

Bx.T = \{Tdl- Px d Px ) 

The bulk dynamics conserve two quantities: the to- 
tal energy H N — XL=i e ' x an< ^ ^fie total deformation 
1Z N = ~^2 x= i t x . The energy conservation law can be 
read locally as ([|, [|) 

e x (t)-e x (0) = J e x (t) - Jx+i(t) 

where J x (t) is the total energy current between x — 1 and 
x up to time t. This can be written as 

J e x (t) = N 2 f f x (s)ds + M x {t) 
Jo 

In the above, M x (t) is a martingale, i.e. a stochastic 
noise with mean 0. The instantaneous energy current j x 
can be written as 

fx = -r x -\Px ~ |v(p 2 ) 

The first term —r x ^\p x is the Hamiltonian contribution 
to the energy current while the noise contribution is given 
by the discrete gradient -(7/2)V(p 2 ) = ("f/2)(p 2 x -p 2 x+1 ). 

Similarly, the deformation instantaneous current j x be- 
tween x — 1 and x is given by 
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~Px 



We denote by fi ss —< ■ > ss the invariant probability 
measure for the process. In the case Tg = T r = T, the 
system is in thermal equilibrium. There is no heat flux 
and the Gibbs invariant measure (or canonical measure) 
is a product Gaussian measure /i ss = fi T ' e depending on 
the temperature T and the mean deformation I: 



N-i 




(2) 
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III. FLUCTUATION-DISSIPATION EQUATION 
AND HYDRODYNAMIC LIMIT 

Diffusive interacting particle systems can be classified 
in two categories: gradient systems and non-gradient sys- 
tems ((HI)- For the first, we can write the curent of the 
conserved quantities as a spatial discrete gradient. For 
example SSEP and KMP process are gradient s yst ems. 
A powerful approach introduced by Varadhan ([16j]) to 
study non-gradient systems is to obtain a fluctuation- 
dissipation equation, meaning a decomposition of the 
current j of conserved quantities as the sum of a mi- 
croscopic gradient V/i and of a fluctuating term of the 
form Cu: 

j = X7h + £u (3) 

where C is the generator of the interacting particle sys- 
tem. In fact, the equality © is only an approximation 
in a suitable Hilbert space f[14|). 

Fortunately, for our system, we can write an equal- 
ity like ((3]) without approximations. The fluctuation- 
dissipation equation for the deformation current j r and 
the energy current j e is given by (0j) 

\3t = V [</> x ] + Cg x [ ! 

where 

1 , 7 , 1 7 „ 

^ = 2^* + 2 Px + Qq >xPx+1 + 4 V( ^+i) 

and 

h -i Vl Vx s 

h x =7 Px, 9x = + 2~V X + r *-V 

Assume that initially the system starts from a local 
equilibrium < • > with macroscopic deformation profile 
uo(q) and energy profile eo(q), q € [0,1]. This means 
that if the macroscopic point q £ [0, 1] is related to the 
microscopic point x by q = x/N then at time t = 

< r[ N q\(G) >^ u Q (q), < e[jv,](0) >^ e (q) 

as N goes to infinity. The currents are related to con- 
served quantities by the conservation law 

dt < r [Nq] {t) >« -Nd q < f [Nq] (t) >, 
dt < e [Nq] (t) >« -Nd q < j( Nq] (t) > . 

By (j3J) and the fact that the terms N < Ch x > and 
N < Cg x > are of order 0(N^ 1 ) and do not contribute 
to the limit (0) we get 

dt < r [Nq] (t) >» 7 -!A < r [Nq] (t) > 
d t < e [Nq] (t) >w A < <f> [Nq ](t) > 



To close the hydrodynamic equations, one has to replace 
the term < </>[jv 9 ] (t) > by a function of the conserved 
quantities < rm q i(t) > and < erjyJt) >. 

The replacement is obtained through a "thermal lo- 
cal equilibrium" statement (see @,@,[I3,[ll,[ll,[l|]) 
for a rigorous justification in the context of conservative 
interacting particle systems). We repeat here the argu- 
ments of [J] for convenience of the reader. Thermal lo- 
cal equilibrium assumption corresponds to assume that 
each given macroscopic region of the system is in equilib- 
rium, but different regions may be in different equilibrium 
states, corresponding to different values of the parame- 
ters. Let us consider an atom with position q = x/N 
which is far from the boundary and introduce a very large 
number 2L + 1 of atoms in microscopic units [L 1), 
but still an infinitesimal number at the macroscopic level 
((2L + l)/N < 1). We choose hence L = eN where 
e <C 1 in order to have these two conditions. We consider 
the system in the box A^ (x) composed of the atoms la- 
beled by x — L, . . . , x + L. The time evolution of the 
2L + 1 atoms is essentially given by the bulk dynam- 
ics; since the variations of deformation and energy in the 
volume containing the 2L + 1 atoms changes only via 
boundary effects and we are looking at what happened 
after TV 2 microscopic time units, the system composed 
of the L atoms has relaxed to the micro-canonical state 
Xf q (t),e (t) corresponding to the local empirical deforma- 
tion f q (t) and the local empirical energy e q (t) in the box 
Al(x). This means that we can divide the observables 
into two classes, according to their relaxation times: the 
fast observables, which relax to equilibrium values on a 
time scale much shorter than t and will not have any 
effect on the hydrodynamical scales and the slow observ- 
ables which are locally conserved by the dynamics and 
need much longer times to relax. We can then replace 
the term < <j>[ Nq ](t) > by *r,(t),e t (t)(<fo)- % equivalence 
of ensembles, in the thermodynamic limit N — * oo and 
then e — > 0, this last quantity is equivalent to 

7+ 2 7 < e [Nq] (t) > + 7 A 1 (< r [Nq] (t) >f 

We have obtained the time evolution of the defor- 
mation/energy profiles u(t,q) = lim < rm q i(t) >, 
e(t,q) = lim < erjvg](i) > in the bulk. At the bound- 
aries, Langevin baths fix temperature at Tg and T r . 
Hence it is more natural to introduce the couple of 
deformation/temperature profiles rather than deforma- 
tion/energy profiles. The temperature profile T{t, q) 
is related to u(t, q) and e(t, q) by e(t, q) — T(t, q) + 
u(t,q) 2 /2. Deformation and temperature profiles evolve 
according to the following equations 

d t T= i(7 + 7- 1 )AT + 7 - 1 (V U ) 2 , 

dtu = 7 -1 Aw, 
' T(t,0)=T e , T(t,l)=T r , ( 5 ) 
u(t,0) = u{t,l) =£, 
T(0,q) = T Q (q), u(0, q) = U (q). 
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As t goes to infinity, the system reaches its steady 
state characterized in the thermodynamic limit by a lin- 
ear temperature profile T(q) = Ti + (T r — Te)q and a con- 
stant deformation profile f(q) — I. The system satisfies 
Fourier's law and the conductivity is given by (7+7~ 1 )/2 

(!)• 



IV. INTERPRETATION OF THE 
FLUCTUATION-DISSIPATION EQUATION 

We have seen that functions h x and g x had no influence 
on the form of the hydrodynamic equations. This is well 
understood by the fact that they are related to first order 
corrections to local equilibrium as we explain below. 

Assume T e(r) = T± ST/2 with ST small. For ST = 0, 
the stationary state < • > ss equals the Gibbs measure fdp 
(see [2]). If ST is small, it is suggestive to try an ansatz 
for < • > ss in the form: 



M = z 1 II d Vx dr * ex P (~ 2T{x/N)^ x + <yVx ~ ^ 



where T(-) is the linear interpolation on [0,1] between 
Ti and T r . \x is the "local equilibrium" approximation of 

< • > ss . Let f ss be the density of the stationary state 

< • > ss with respect to fl, i.e. the solution of C*' T ^ 'f ss — 
0. Here £*' T (') is the adjoint operator of C in L 2 (/i). It 
turns out that 



£*,T(.) = _ A + lS + Bl , 



D 



N,T r 



N-2 



N-l 



T2 I M / j -Jx,x+1 *- / j Jx,x+l 



x=l 
N-l 



ST. 



Prl>, ■ i A.,.,.. : I + —(dp, - d'p N ) 



x=l 

\2\ , mi *j-l\ 



+ 0{{STY) + 0(N- 

where j e and j r are the energy and deformation currents 
for the reversed dynamics at equilibrium. They are ob- 
tained from j e and f by reversing momenta p — > —p. 
Expanding / ss at first order f ss — 1 + ST v + o(ST) , we 
get that for large N and small gradient temperature ST, 
v has to satisfy the following Poisson equation: 



{-A + ~fS)v =T~ 



N 



N-2 
x = l 



N-l 



i 

N 



Let v the function obtained from v by reversing momenta. 
By the fluctuation-dissipation equation (j4|) we get 



N-l 



NT 2 



^(gx-^ + OiN- 1 ) 



Therefore the functions g x and h x are directly related to 
first order corrections to local equilibrium. 



V. NON-EQUILIBRIUM FLUCTUATIONS AND 
STEADY STATE CORRELATIONS 



Assume that initially the system starts from a local 
equilibrium < • > with macroscopic deformation profile 
uo(q) and temperature profile T (q), q G [0,1]. The time- 
dependant deformation fluctuation field and energy 
fluctuation field Y t N are defined by 



1 



iV 



V X— 1 



where H, G are smooth test functions, (T(t, ■), u(t, •)) are 
solutions of the hydrodynamic equations ([5]) with e = 
T + u 2 /2. 

The fluctuation-dissipation equations (section llllj) give 
(0, H): 



R?(H) = R»(H) + -J* R?(AH)ds + M 



UN 



Y t N (G)=Y N (G)+ 1 J Q t Y s N (AG)d 



So ds 



+ Mt 



2.N 



AT 



Exer N (&G)(x/N)f x (u 



where ftA 1,N and A4 2 ' N are martingales and f x is the 
function defined by 



7 



V 2 
-T 



— Px + lPx - -f P 3 



Covariance of the limit martingales are computed us- 
ing standard stochastic calculus and thermal equilibrium 
property (fl, 0): 



A4 



UN 



2 /■* 
- / ds 

7 Jo 



[0,1] 



dqT(q,s)(VH) 2 (q) 



M 



2,N 



dq / d™ 2 ( g , s )T( 9 , s )(VG) 2 ( 9 ) 
7J[o,i] Jo 

+( 7 + 7 _1 ) ! ds f dqT 2 (q,s)(VG) 2 (q) 
Jo -'[0,1] 



Ml' N M 2 t > N 

- [ ds [ dqu(s,q)T(s,q)(WG)(q)(WH)(q) 
7 Jo J[o,i] 

Hence R^ converges as N goes to infinity to the solu- 
tion of the linear stochastic differential equation: 

1 



d t R = -AR - V 

7 



-T(t,q)W 1 (t,q) 



7 



(6) 
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where W\{t, q) is a standard space time white noise. 

The description of the limit for the energy fluctuation 
field is more demanding. We have first to close the equa- 
tion. In order to do it, we use a "dynamical Boltzmann- 
Gibbs lemma" ([HI], [H). Observables are divided into 
two classes: non-hydrodynamical and hydrodynamical. 
The first one are non conserved quantities and fluctu- 
ate in a much faster scale than the others (in the time 
scale where these last change). Hence, they should aver- 
age out and only their projection on the hydrodynamical 
variables should persist in the limit. One expects there 
exist constants C, D such that 



N 



d S ^(AG)(x/7V){/ x (o 



x = l 



-C{r x - u(s, x/N)) -D(£ x - e(s, x/N))} 



vanishes as N goes to infinity. Constants C and D de- 
pend on the macroscopic point q = x/N and on the time 
t. In order to compute these constants, we assume ther- 
mal local equilibrium. Around the macroscopic point q, 
the system is considered in equilibrium with a fixed value 
of the deformation u(t, q) and of the temperature T(t, q). 
The constant C, D are then computed by projecting the 
function f x on the deformation and energy fields ([01], 
If fi T ' e is the Gibbs equilibrium measure with tem- 
perature T and mean deformation £ (the mean energy is 
then e = £ 2 /2 + T), we have $(l,e) = ^(/ x ) = e+£ 2 /2 
and then 

C = di<S>(u(s, q), e(s, q)), D = d E <S>(u(s, q), e(s, q)) 

Therefore the time-dependant energy fluctuation field 
Y t converges as N goes to infinity to the solution of the 
linear stochastic differential equation: 



d t Y 









( 7 + i) Ay h 


■i 


V7 7 ) 









^J 1 + 1 -^T{t, q)W 2 {t, q) + u(t, qJ^l^-W^t, q) 



(7) 



r 



where W2(t, q) is a standard space-time white noise in- 
dependent of W\{t, q). 

Remark that the deterministic terms in ((6]) and ([7|) 
result from linearizing the nonlinear equation as ([5]). 

We now compute the fluctuations fields for the NESS 
< • > ss which is obtained as the stationary solution of 
the Langevin equations ([B][7]). The field L t defined by 
L t = —£Rt + Y t is solution of the Langevin equation 



d t L = bAL - V 



2bT(t,q)W 2 (q,t) 



1 



with b = -(7 + 7 1 ). The fields Rt and L t are solutions 

of independent decoupled linear Langevin equations and 
converge as t goes to infinity to independent Gaussian 
fields. It follows that Rt and L t converge to stationary 
fluctuation fields R ss and Y ss such that 



Cov{R ss {G),R ss (H))= / dqG(q)H(q)T(q) 



Cov(Y ss (G),Y ss (H)) = / dqG{q)H{q){f 2 {q)+£ 2 f{q)} 
Jo 

+2(T i ~T r ) 2 f G(q)(A- x H){q)dq 
Jo 

Cov(Y ss (G),R ss (H))=t f H(q)G(q)f(q)dq 

Jo 

Observe that the covariance of the fluctuations of en- 
ergy is composed of two terms. The first one corresponds 



to Gaussian fluctuations for the energy under local equi- 
librium state while the second term represents the con- 
tribution to the covariance due to the long range corre- 
lations in the NESS. As in the case of SSEP and KMP 
process, the correction is given by the Green function of 
the Dirichlet Laplacian (0], (20|). 



VI. 



LARGE FLUCTUATIONS 



A. Macroscopic dynamical behavior 

Assume that initially the system is prepared in a state 
with a deformation profile uo, energy profile £0 and hence 
temperature profile Tq = Eq — u^/2. In a diffusive scale 
the deformation (resp. energy, resp. temperature) pro- 
files u (resp. e, resp. T) where e = T + u 2 /2 evolve 
according to the hydrodynamic equations ([5]). 

Our aim is to obtain the large deviation principle 
corresponding to the law of large numbers ([5]) . It consists 
to estimate the probability that the empirical quantities 
(deformation, energy, temperature) do not follow the 
corresponding solutions of (O but remain close to some 
prescribed paths. This probability will be exponentially 
small in N and we look for the exponential rate. We 
follow the classic procedure in large deviation theory 
(Q,[l3|): we perturb the dynamics in such a way that 
the prescribed paths become typical and we compute 
the cost of such perturbation. 
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Fix a path y(t, 
formation profile and empirical energy profile S,f are 



(u(t, -),e(t, •)). The empirical de- 

cN 
"t 



defined by 

11? (q) = N- 1 J2 r x (t)l [x/N< ( x +x)/N)(q), (8) 

x=l 
N 



N 



In appendix, we explain how to define a Markovian dy- 
namics associated to a couple of functions H(t, q),G(t, q), 
q G [0, 1], such that the perturbed system has hydrody- 
namic limits given by u and e. This is possible if the 
function F = (H, G) solves the Poisson equation 



d t y = Ay-V(aVF)) 
F(t,0) = F(t,l) = (0,0) 



(9) 



where the mobility a := a[u,e) is given by 

a(u, e) = 2 ( ? T u2 ^ T + T2 V T = E - u 2 /2 (10) 

The perturbed process defined a probability measure 
P on the deformation/energy paths space by mean of the 
empirical deformation and energy profiles (see ©). 

Our goal is to estimate the probability 

,N cN\ 



= E 



F[(R.?,£?)~(u(s r ),e(a,-)), sG[0,i\] 



{(■Rf ,£f )~(r(«,-),e(«,0), »e[0,t]} 



To avoid irrelevant complications due to the fluctua- 
tions of the initial state which have no incidence on the 
derivation of the quasi-potential, we assume that the ini- 
tial profiles uq and To are the stationary profiles f(q) = £ 
and f(q) =T e + (T r - T t )q. The function F is such that 



(«(«.-), »e[<M]] 



l 



In the appendix we show that in the large N limit, under 
P, the Radon-Nikodym derivative is given by 



w exp {-NJ[ ,t] [u, e) j 



where 



1 f* 

J[0,t]( u > £ ) = 2j ds<WF( S ,-),aWF( S ,-)> q (11) 

where a is here for a(u(s, •), e(s, •)) and < •, • > q for the 
usual scalar product in L 2 ([0, 1], dq). Hence we have ob- 
tained 

P[(ft^£jV(u(v), £(*,•)), *€[0,t]] 
w exp {-7VJ[o !t ](w,e)} 



B. The quasi-potential 

To understand what is the quasi-potential, consider 
the following situation. Assume the system is macro- 
scopically in the stationary profile (u(— oo, •), e(— oo, •)) = 
(I, f (•) + £ 2 / 2 ) at * = -oo but at f = we find it in the 
state (u(q),e(q)). We want to determine the most prob- 
able trajectory followed in the spontaneous creation of 
this fluctuation. According to the precedent subsection 
this trajectory is the one that minimizes J[_oo,o] among 
all trajectories (u, e) connecting the stationary profiles to 
(u,e). The quasi-potential is then defined by 



W(u,e) 



inf J[ 0j t](u,e) 

(u,e) 



MFT postulates the quasi-potential W is the appropriate 
generalization of the free energy for non-equilibrium sys- 
tems and this has been proven rigorously for SSEP ([3(). 
W is solution of an infinite-dimensional Hamilton- Jacobi 
equation which is in general very difficult to solve. It has 
been solved for specific models (SSEP and KMP) having 
a single conservation law ( Q , [J] ) . For the system we con- 
sider, two quantities are conserved and we are not able to 
solve this Hamilton- Jacobi equation. Nevertheless we can 
compute the quasi-potential for the temperature profile 
(TrZ|) in the case 7 = 1. The latter is obtained by pro- 
jecting the quasi-potential W on the deformation/energy 
profiles with a prescribed temperature profile. 

Consider the system in its steady state < • > ss . Our 
aim is here to estimate the probability that the empirical 
kinetic energy defined by 



S N (q) = N- 



N 

1 ^2pl 1 [x/N,(x+l)/N)(q) 



(12) 



is close to some prescribed temperature profile 7r(q) dif- 
ferent form the linear profile f(q) = T e + (T r - T t )q. This 
probability will be exponentially small in TV 

([e%)~7r(g)]) ss «exp(-W(7r)) 

By MFT, the rate function V(w) coincides with the fol- 
lowing projected quasi-potential 

V{w) = M inf J [0tt] (u,e) 

where the paths set At^ is defined by 

u 2 (t,-) 



<At.7r 



(u,e); e(f,-) 



»(o,-) = 4 r(o,-) = T(.) 



Paths y — (u, e) <G At.n must also satisfy the boundary 
conditions 

u(t, 0) = u(t, 1) = £, e(t, 0) = T e +£ 2 /2, e(t, 1) = T r +£ 2 /2 

(13) 
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In fact, it can be shown that Jm(\(u,e) — +00 if the 
path y does not satisfy these boundary conditions. 

Our main result is the computation of the projected 
quasi-potential: 



V(tt) = inf^TT, r) 



(14) 



where T = {r e C x ([0, 1]); r'(q) > 0, r(0) = T t , r(l) = 
T,.} and 



•H 71 ", r ) = / d 1 
Jo 



— - 1 - log — -r^- - log 



t(«) b (T r - 2i) 



Before proving (|14[) let us make some remarks. First, 
VJw) is equal to the rate function for the KMP process 
(Q). Nevertheless, it is not easy to understand the deep 
reason. The symmetric part S of the generator C is more 
or less a time-continuous version of the KMP process for 
the kinetic energy but the Hamiltonian part has a non- 
trivial effect on the latter since it mixes momenta with 
positions. Hence, the derivation of the quasi-potential 
for the kinetic energy can not be derived from the the 
computations for the KMP process. Secondly, we are 
able to compute V only for 7 = 1. When 7 is equal to 
1 hydro dynamic equations for the deformation and for 
the energy are decoupled but since temperature is a non- 
linear function of deformation and energy, it is not clear 
why it helps- but it does. Formula ([T3]l shows that the 
large deviation functional V is nonlocal and consequently 



not additive: the probability of temperature profile in 
disjoint macroscopic regions is not given by the product 
of the separate probabilities. Nonlocality is a generic 
feature of NESS and is related to the 0(N^ 1 ) corrections 
to local thermal equilibrium. 

Let us call S(ir) the right hand side of equality (fT4|) . 

For every time independent deformation/energy pro- 
files (r(q), e(q)) and r(q) £ T we define the functional 

f 1 (T T t' {r-lf 

U(r, e, t) = / dq i 1 - log - - log — + — 

Jo [t t T e -T r 2t 

(15) 

where T(q) — e(q) — r(q) 2 /2 the temperature profile 
corresponding to (r(q),e(q)). Define the function r := 
r(r, e) of T as the unique increasing solution of: 



At 1 
T 2 = t-T- -(r- £) 2 

(Vr) 2 2 [ ' 

r(Q)=T t , r(l)=T r 



(16) 



Fix deformation/energy paths satisfying boundary con- 
ditions (flU|) and define Z by 



y=["), Z=[d t y-Ay + V(aV(6U))}, (17) 



In the appendix we show the following formula 



J [0lt ](«, e) = U(u(t, ■), e(t, •), r(e{t, -),u(t, ■))) - U(u(0, •), e(0, -),t{u(0, ■), e(0, •))) (18) 



where 

Sr 

5U = (u,e,t(u,e)) 
5U 

If (u, e) belongs to -Ai i7r , 

Cr( tt (O,0,e(O,-),r(u(O,-),e(O,-))) 
= U(t,T + £ 2 /2 7 t(£, T + f/2) ) = 

and 

U{u{t, •), e(t, •), r(u(t, •), e(*. •))) > ■^( 7r . ■).«(*. '))• 

The two last terms on the right hand side of (fT8f are 
positive so that for every paths in At.n, we have 

J[o,t](«,e) > S'l 71 ") 



and we obtain hence 



V(ir) > S(ir) 



(19) 



To obtain the other sense of the inequality, we have to 
construct an optimal path (u*,e*) S At,n such that the 
two last terms in the right hand side of (fT5| are equal to 
0, i.e. 



d t y = Ay - V(o-VOH/)) 

u{t,q)=l 



(20) 



We note T* = e* — u* 2 /2 the corresponding temperature. 
Then reporting in (fT5)) . we obtain 

J [0 ,t](ti*,e*) = I7(tt*(i),e*(i),r(u*(t),e*(i))) (21) 
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By the definition (|15|) of U and by using the fact that 
u*(t, q) — £, we obtain 

J m (u*,e*)=HT*(t,-),T(u*(t,-),e*(t,-)) 

= F(tt,t(£,tt + £ 2 /2)) (22) 

The variationalproblem defining S is solved for r = 
t(£,tt + £ 2 /2) (0|) so that 

5(tt) =^(7r,r(£,7r + ^ 2 /2)) 

and therefore we have 

V(tt) = inf inf J [0lt ](«,e) < S(tt) 

This inequality with ((TSJ| shows that V(7r) = S^tt). ft 
remains to prove that such "good" path exists. The proof 
is similar to [H and we shall merely outline it. Equation 
([2H)l is equivalent to the following one 



d t T* = — A(T*) 
u*(t, ? ) =1 



2V 



*\2 



(T*) 



■V(r*) 



(23) 



where r*{t,-) = r{£,T*(t,-) + £ 2 /2). Let us denote by 
0*( S) •) = T*(* - s, •) the time reversed path of T*. 9* 
can be constructed in the following procedure. We define 
0*(s,q), ee[0,t], qe [0,1] by 



'( s ,.) = p( s ,.)-2p( S ,.) 2 



where p(s, q) is the solution of 



Ap(s, 



[(Vp)( S ,-)] 2 



p(s,0)=T t , p(a,l)=T r 
^p(0,q)=p (q)=T(£,7r + e 2 /2)(q) 

It can be checked that T*(s, q) = 9*{t- s, q) solves (|23|) . 
Moreover, we have T*(0, •) = 9*(t, •) and T*(t, •) = tt(-). 
This path belongs to At.* only as t — > oo since 6?*(i, •) 
goes to ?(•) as f -> +oo. We have hence in fact to 
define T* by the preceding procedure in some time in- 
terval [ti,t] and to interpolate T(-) to n*(ti,-) in the 
time interval [0,ti] (see Q, [1] for details). This opti- 
mal path is also obtained as the time reversed solution 
of the hydrodynamic equation corresponding to the pro- 
cess with generator C*. It is easy to show that this last 
hydrodynamic equation is in fact the same as the hy- 
drodynamic equation corresponding to C. This is the 
"generalized" Onsager-Machlup theory developed in Q 
for NESS: "the spontaneous emergence of a macroscopic 
fluctuation takes place most likely following a trajectory 
which can be characterized in terms of the time reversed 
process." Observe also the following a priori non trivial 
fact: the optimal path is obtained with a constant defor- 
mation profile. 



VII. CONCLUSIONS 

In the present work we obtained hydrodynamic limits, 
Gaussian fluctuations and (partially) large fluctuations 
for a model of harmonic oscillators perturbed by a con- 
servative noise. Up to now MFT has been restricted to 
gradient systems with a single conservation law. This 
work is hence the first one where MFT is applied for a 
non-gradient model with two conserved quantities. The 
quasi-potential for the temperature has been computed 
in the case 7 = 1 and it turns out that it coincides with 
the one of the KMP process. Our results show this sys- 
tem exhibits generic features of non-equilibrium models 
: long range correlations and non-locality of the quasi- 
potential. 

Nevertheless our study is not completely satisfactory. 
It would be interesting to extend the previous results to 
the case 7^1 and to compute the quasi-potential for the 
two conserved quantities and not only for the tempera- 
ture. The difficulty is that there does not exist general 
strategy to solve the corresponding infinite-dimensional 
Hamilton- Jacobi equation. 



APPENDIX A: THE PERTURBED SYSTEM 

Let us denote by u)(s) — (p x (s),r x (s)) x the configu- 
ration of the process at time s and P the probability 
measure on the deformation/energy paths up to time t 
that the process (w(s))o< s <t defines. It is well known 
([HI) that if M is a P- martingale with quadratic varia- 
tion [M] t then the probability measure P with Radon- 
Nykodim derivative given by 

^ = exp(M 4 --[M] ( j (Al) 

defines a Markov process (&(s))o< s <t- I n particular for 
a time dependant function f(t,u>) on the configuration 
space, the process 

M t = f(t, w(f))-/(0, w(0))- / ds(d s +N 2 C)f {a, io(s))ds 



is a P-martingale with quadratic variation [M] t given by 

[M] t = N 2 /V(/./)](w(*))ds (A2) 
Jo 

where the "carre du champ" operator •) is defined by 
r(u, v) = L(uv) — uCv ~ vCu 

N-l 

= X x>x+ i(u)X x>x+ i(v) 

x=l 

+ T e (d Pl u)(d Pl v) + T r (d PN u)(d PN v) 

The Markov process Q has then generator given by 

C = C + T(f,-) (A3) 

The first order partial differential operator T(/, •) can 
be seen as a perturbative drift. 
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1. Hydrodynamic limit 

Fix a couple of smooth functions H(t,q),G(t,q) van- 
ishing at the boundaries and define 



By Taylor expansions, one has, up to 0(N 2 ) correc- 
tions, 



/ = fH + fa 

JV-2 N-2 



x=2 



x=2 



where h x , g x are the functions appearing in the 
fluctuation-dissipation equation (f5|). 
By taylor expansion, we get 



X g , z +i(f H ) = -^(VH)(t,z/N)(p z+1 -p z )+0(N- 2 ) 
and 



*,,,+i(/g) = -^(VG)(t, z/N)6 z + o(n- 2 ) 



where 9 Z is defined by 



T(/,M « --(VH)(t,x/N){p x+1 (p x+1 -p x ) 

-Px-l(Px -Px-l)} 



-(VG)(i, x/N) { Px+l 9 x - p x -i9 x ^} 



and 



r(/, 9x) « -— (Vff s/JVT) {(r* + r,_i) 

(Px+l + Px-1 -PxPs-l)} 



27V 



{VG)(t,x/N) {{r x + r x -x){B x p x+l - ^-ipx-i)} 



#2 = PzPz + l + ^(Pz+l^z + 

By (UJ) and (|A3|) . the instantaneous deformation (resp. 
energy) current x (resp. B ) for the perturbed 

system is now 

[f x = -V{r x ) + Lh x -T{f, h x ) 

| (resp. f x = Vcj) x +Cg x -T(f,g x ) - p x p x -iX x _ hx {f) ) 

In comparison with fluctuation-dissipation equation |4|. 
currents are modified by terms of order 0(N~ 1 ) (see be- 
low). 



and 



p x -ip x X x - ljX f « -— (VH)(t,x/N)p x ^ip x (p x -p x -i) 



(VG)(t,x/N)d x ^ 1 p x _ lPx 



N 



To obtain hydrodynamic equations for the perturbed 
process Cj we use thermal equilibrium property. Observe 
that at equilibrium under the Gibbs measure /i^ with 
mean deformation I and temperature T, we have 



^ T " l {^{f,9x) +p x -ip x X x „ lyX f) 



-N- l a(tT + e/2)( ^ H )^ X I N ) 



with the mobility matrix cr defined in (| 10[) . We repeat 
the arguments of section Mil and we get that the hydro- 
dynamic limit of Cj is given by 

d t y = Ay- v(aVF)) 

with boundary conditions like in 

2. The dynamical large deviations function 

Fix a path 3^ = (u, e) and consider the perturbed pro- 
cess defined above with F = (H, G) chosen according to 



([9]). We show here that in the large N limit 



E 



dP 
-^1 



{(K? ,£f )~(«(a,-),e(a,-)), a6[0,t]} 



with J[ 0it] (u,e) defined in ([TT)). By (|A"l~j) and ijMj) . 
the logarithm j7[o,t]( ti; ) of the Radon- Nykodim derivative 
dP/dP is given by 



Jio, t] M = /(i,w(t))-/(0,w(0)) 
TV 2 /■* 

d s (d s + N 2 £)f( Sl La( S ))d S + — [T(fJ)](Lo(s))ds 
o z Jo 
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Observe that 

C(r x + Vh x ) =Ar a + V(r(/,/i x )) 

and 

L(e x + Vg x ) = A<f> x + V(r(/, g x ) + p x -ip x X x _i, x f) 

Moreover the term T(f, f) = T(f, f H ) + T(f, f G ) gives 
a contribution equal to 

N~2 



T{fJ)=Y J H{t,x/N)\/T{h x ,f) 

x=2 

J2 G(t, x/N)V [T{g x , /) + p x -ip x X x _ liX f] 



N-2 



Recall that F = (H, G) has been chosen such that 

p[{(7ef,£f)~( u ( v ), £(«,•)), se[o,t}}] =1 

By local equilibrium statement, integration by parts and 
the precedent computations for T(f,h x ),T(f,g x ) and 

Px — lPx X x — 1 x /, one has 

N f* 



so that 



E 



3^1{(7Jf ,£f )~(«(s,-),e(s,-)). s£[0,t]} 



= -iVJ [ot] («,£) 



APPENDIX B: PROOF OF FORMULA (fl8]> 

The goal is to express J[ .t](u,e) as the 
sum of U(u(t,-)),e(t,-),T(e(t, -),u(t,-))) 
U(u(0, -)),£(0, •),r(e(0, -),u(0, •))) and positive terms. 
We recall that 

1 /■* 

Jm(«,£) = 2^ ^(Vf( S ,.) lff VF(v)}, (Bl) 

where a = a(u(s, •), s(s, •)) is the mobility matrix defined 
in pop . By the definitions of y and Z (see (|17p). we have 

V[crVF] = Ay - 9t^ = V(ctV(<5[/)) - Z 

Inserting this last expression in (|B1[) we get 



J m (u,e) = ijf d» < V- 1 [V(<rV((St7)) - Z] . 

cr^V" 1 [V(crV((5t/)) - Z] > g 

We develop the expression and we get 

J m (u,e) = X -l ds(V- 1 Z,a- 1 V- 1 Z) q 
1 I* 

+ - ds< aV(SU),V(SU) > q 
* Jo 

~ f ds<V~ 1 Z,V(5U)> q 
Jo 



By integration by parts and because of the boundary 
conditions, the second term and third term on the right 
hand side are equal to 

[ ds< d s y, SU > q ~ f ds< aV(SU),V(6U) > q 
Jo 2 J 

- f ds< Ay, SU) > g 
Jo 

The first term is the integral of the time derivative of 
s — > U(u(s, ■), e(s, •), t(u(s, ■), e(s, •)) because 



SU, 



(r,e,r(r,e)) = 



Hence it is equal to U(u(t, •), e(t, •), r(u(t, •), e(t, •)) — 

?7(u(O,0,e(O,-),T(tt(O,-),e(O,-)). 

We develop the two other terms using the expression 
of SU. A simple computation shows 



SU 



Sr 

Hence we get 



r/T - £/t, 



SU 



T -i 



1 r* 

- ds< aV(SU),V(SU) > 

2 Jo 



, (Vu) 2 (VT) 2 VtVT _(Vt) 
o J [o,i] \ T T 2 t 2 t 



VuVt 

+2(1 -u) 5— + T(u-l)- 



,(Vr 



For the term 



/ ds<Ay,6U) > q 
Jo 



we perform an integration by parts and we obtain 



1 f f(Vu) 2 VuVt 
ds dq{ + £- 

o J [0,1] 



T 



(B2) 



(B3) 



(B4) 



(VT) 2 VTVt VuVt 



u- 



T 2 T 2 T* 

The sum of the expressions (|B2p and (|B4|) is equal to 



VrVT , (Vr) 2 
— + T 2 — J- 



ds / dq 

o J [o,l] I T 

VuVt 

+ (l-u) — + T{£-u)- 



Remark now that by integration by parts, 



/ dq\7(T-r)^ + f 
J [0.1] T J\l 



, (Vr) 2 
dq 5- 

[04] t2 



VtVT 
dq 2 — 

[0,1] T J [0,1] 

[0,1] \ T J J [OA] 



/ dq(T-T)(- 
J[o,i] V 
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and 

f , , VuVt If n9i Vt 
/ dq(i-u)—z-=-~ dqV[(£-u) 2 — 
J lo,i] ' 2 J[o,i] T 

1 f , ,JAt (Vt) 2 1 

Collecting all these facts and using the equation defining 
r, we obtain (fl~8)) . 
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